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Abstract 

Dilute Bose gases, cooled down to low temperatures below the Bose-Einstein con- 
densation temperature, form coherent ensembles described by the Gross-Pitaevskii 
equation. Stationary solutions to the latter are topological coherent modes. The 
ground state, corresponding to the lowest energy level, defines the standard Bose- 
Einstein condensate, while the states with higher energy levels represent nonground- 
state condensates. The higher modes can be generated by alternating fields, whose 
frequencies are in resonance with the associated transition frequencies. The conden- 
sate with topological coherent modes exhibits a variety of nontrivial effects. Here it 
is demonstrated that the dynamical transition between the mode-locked and mode- 
unlocked regimes is accompanied by noticeable changes in the evolutional entangle- 
ment production. 
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1 Introduction 

Bose-Einstein condensation, as is known, occurs when a macroscopic number of bosons 
piles down to the same lowest single-particle state. In equilibrium, this is the sole possible 
type of Bose-Einstein condensate. When the Bose-condensed system is rotated, vortices 
can appear. The Bose-Einstein condensate with vortices is the most known example of a 
nonground-state condensate. At low temperatures Bose gases are described by the Gross- 
Pitaevskii equation (see reviews [1-3]). In the presence of a trapping potential, there arises 
a whole spectrum of energy levels associated with the stationary solutions to the Gross- 
Pitaevskii equation. Each of these stationary solutions is, by definition, a topological 
coherent mode [4]. The lowest energy level corresponds to the standard Bose-Einstein 
condensate, while the higher states describe various nonground-state condensates. The 
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latter can be generated by alternating fields, whose frequencies should be in resonance with 
the related transition frequencies [4-8]. It is feasible to realize an oscillatory modulation 
of either the trapping potential or of the atomic scattering length. The trapping potential 
can be either single-well or multiwell. 

Nonground-state Bose-Einstein condensates, with generated topological coherent 
modes, exhibit several interesting effects, such as interference fringes, interference current, 
mode locking, dynamical transition, critical phenomena, chaotic motion, harmonic genera- 
tion, parametric conversion, and atomic squeezing [4-8] . The aim of the present communi- 
cation is to study the behaviour of entanglement production in a trapped nonground-state 
condensate with topological coherent modes and to demonstrate that the properties of this 
entanglement production experience noticeable changes under the dynamical transition 
from the mode-locked regime to mode-unlocked regime. 



2 Topological Coherent Modes 

Dilute Bose-condensed gases at low temperature are described [1-3] by a coherent field 
<p(r,t) satisfying the Gross-Pitaevskii equation 

ih^<p(T,t) = (H[p]+V) <p(T,t), (i) 

in which V = V(r,t) is a modulating potential and the nonlinear Hamiltonian 

H[<p] = -h 2 ^ + U(r) +NA S \<p(r, t)\ 2 (2) 

contains a trapping potential U(r) and the interaction intensity A s = AiiTi 2 a s /m with m 
being atomic mass; a s , scattering length; and N, the total number of atoms. 

The topological coherent modes are the solutions to the stationary Gross-Pitaevskii 
equation 

H[ip n ] ip n (r) = E n ip n (r) , (3) 

where n is a labelling multi-index. The generation of topological modes is accomplished 
by transferring atoms from an occupied mode (say, ground state) to a higher mode by 
means of an alternating field 

V(r, t) = - B{r)e iujt + - B*(r)e~ lu)t , (4) 

whose frequency uj is in resonance with a chosen transition frequency u>2i = (E 2 — Ei)/K 
between two selected energy levels, implying that the detuning Aa; = uj — ou^i is small, 
|Aa>| uj. Generally, it is admissible to apply several resonant fields of form (4) con- 
necting different modes. The procedure of such a resonant mode generation has been 
thoroughly explained [4-8]. 
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The solution to Eq. (1) is expressible as 

¥>( r »*) = ^c n (t) <p n (r) exp (- ^ E nt^j , (5) 

where c n (t) are functions of time to be found. Substituting expansion (5) into Eq. (1) 
makes it possible to obtain the equations for c n (t). The dynamics of two selected modes, 
connected by the resonant field (4), is described by the temporal dependence of the 
corresponding coefficients C\{t) and Oz{t), whose moduli squared, \c n \ 2 , define the fractional 
mode populations. Let us introduce the population difference s = | C2 1 2 — |ci| 2 and the 
phase difference x = i ln(ci | C2 1 / 1 ci | C2 ) . 

The dynamical properties of a given atomic system, prepared in a state with fixed 
initial conditions, can be governed by varying the amplitude \/3\ = (ipi, Bip 2 )/h and the 
detuning Auj = uo — uj 2 \ of the modulating field (4), with (</?i, ^2) being a scalar product. 
It is convenient to deal with dimensionless parameters normalizing \/3\ and Auj by the 
value a ~ A s N/h characterizing the atomic interaction strength. Then one deals with 
the parametric manifold {b, 5}, in which b = \P\/a is the dimensionless modulating-field 
amplitude and 5 is the dimensionless detuning. On this parametric manifold, there is a 
critical line {b c , S c } separating two qualitatively different regimes of motion. From one 
side of the critical line, there occurs the mode-locked regime, when the fractional mode 
populations are locked in their initial half-planes, such that for all times t one has either 
< |c n (t)| 2 < 1/2 or 1/2 < |c n (t)| 2 < 1, depending on the initial value c n (0). And from 
another side of the critical line, the mode populations become unlocked, oscillating in the 
whole region < \c n (t) 2 < 1, which corresponds to the mode-unlocked regime. In this way, 
by varying either b or 5, and crossing the critical line {b c , 5 C }, one can realize the dynamic 
transition from one regime of motion to another, between the mode-locked and mode- 
unlocked regimes. On the plane of the variables s and x, the dynamic transition happens 
when the initial point of a trajectory crosses a saddle separatrix given by the equation 
s 1 j1 — b\Jl — s 2 cosx + 5s = \b\. Thus in Fig. 1 the initial point {—1, 0}, for the parameter 
b < b c , is below the saddle separatrix, because of which the total related trajectory lies 
below the separatrix, which corresponds to the mode-locked regime of motion. In Fig. 2 
the same initial point, but for b > b c , is already above the separatrix, which results in the 
mode-unlocked regime. 

Figure 1: Phase portrait on the plane{s,x} for a two-mode condensate, with the 
parameters 5 = —0.1, b = 0.35; b c = 0.39821. 

Figure 2: Phase portrait on the plane{s,x} for a two-mode condensate, with the 
parameters 5 = -0.1, b = 0.49; b c = 0.39821. 

3 Evolutional Entanglement Production 

The dynamic transition, with changing regimes of motion, should also be accompanied 
by noticeable changes in all dynamical characteristics of the system, for instance, in 
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the behaviour of atomic squeezing [6]. Here we concentrate on the temporal behaviour 
of entanglement production in the considered system. Entanglement is an important 
characteristic that is assumed to be exploited for quantum information processing and 
quantum computation. We consider here the entanglement generated by the second- 
order density matrix p 2 for the given coherent two-mode condensate. The entangle- 
ment generated by p 2 can be described [9,10] by the entanglement-production measure 
e 2 (t) = e(p 2 ) = log(||p2| Id/Hp® ||d), in which the logarithm is to base 2, pf is the nonen- 
tangling counterpart of p 2 and D is the set of disentangled states formed of the products 
<Pm( r i) < Pn( r 2)- For the density matrix p 2 , we obtain £2(0 = — logsup n |c n (t)| 2 . 

Numerical calculations for the entanglement-productions measure s 2 (t) are presented 
in Figs. 3 and 4 for the initial conditions Cx(0) = 1 and c 2 (0) = 0. Fixing the zero 
detuning, we vary the modulating-field amplitude b. In Fig. 3, the latter is below the 
critical value b c , while in Fig. 4, it is above b c . As is seen, the entanglement production 
is rather different in the mode-locked regime, when b < b c , and in the mode-unlocked 
regime, when b > b c . The considered entanglement, generated in a trap, quantifies the 
amount of correlations between each pair of atoms. If atoms are released from the trap, 
say in the time-of-flight regime, then the studied measure e 2 {t) will change in time, with 
the entanglement produced in a trap serving as the initial condition. Depending on the 
latter, the pair correlations of outcoupled atoms will be different and could be regulated by 
preparing the required entanglement production of trapped atoms. A two-mode system 
is mathematically similar to two-level atoms or to 1/2 spins. Therefore such a system 
can also be employed for realizing quantum gates, which can be done by switching on 
and off the external field. Thus, by varying the amplitude of the modulating field one 
can essentially modify the characteristics of entanglement production, which could be 
employed for information processing. 
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Figure 1: Evolution of the entanglement produced by the second-order density matrix for 
5 = 0, b = 0.3. Time is measured in units of a. 

In conclusion, the nonground-state Bose-Einstein condensates, with resonantly gen- 
erated topological coherent modes, are the objects whose dynamical properties can be 
straightforwardly regulated by external modulating fields. By means of such fields, it is 
feasible to switch the dynamic behaviour between two qualitatively different regimes, the 
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Figure 2: Entanglement production of the second-order density matrix as a function of 
time, measured in units of a, for 5 = 0, b = 0.7. 

mode-locked and mode-unlocked regimes. This change-over from one regime to another al- 
lows one to radically alter the system evolution, with the related changes in the behaviour 
of the fractional mode populations, atomic squeezing, and entanglement production. 
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